Soft-Gluon Resummation and the Valence Parton Distribution Function of the Pion 
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We determine the valence parton distribution function of the pion by performing a new analysis 
of data for the Drell-Yan process tt~N — > /i + fj,~X. Compared to previous analyses, we include 
next-to-leading-logarithmic threshold resummation effects in the calculation of the Drell-Yan cross 
section. As a result of these, we find a considerably softer valence distribution at high momentum 
fractions x than obtained in previous next-to-leading-order analyses, in line with expectations based 
on perturbative-QCD counting rules or Dyson- Schwinger equations. 



PACS numbers: 12.38.Cy,12.38.Bx,14.40.Be 

Although the pion is one of the most important par- 
ticles in strong-interaction physics, our knowledge about 
its internal quark and gluon "partonic" structure is still 
rather poor. Most of the available information comes 
from Drell-Yan dimuon production by charged pions in- 
cident on nuclear fixed targets [3-0| . These data primar- 
ily constrain the valence distribution v n = = d% = 
dl = tZj . Several next-to-leading order (NLO) analy- 
ses of the Drell-Yan data have been performed [4j-|6[. A 
striking feature has been that the resulting valence dis- 
tribution v w (x, Q 2 ) turned out to be rather hard at high 
momentum fraction x, typically showing only a linear 
(1 — x) 1 or slightly faster falloff. This finding is at vari- 
ance with predictions based on perturbative QCD 7J, 
and calculations using Dyson-Schwinger equations 8|, 
for which the falloff is expected to be ~ (1 — x) 2 . On 
the other hand, Nambu-Jona-Lasinio [t| and constituent 
quark models (loj , as well as duality arguments (llj , favor 
a linear behavior. The high- a; behavior of is widely 
regarded to be an important so-far-unresolved problem 
in strong- interaction physics • 

In the kinematic regimes accessed by the fixed-target 
Drell-Yan data, perturbative corrections beyond NLO 
may be significant fl3| ] . The relation z = Q 2 / X1X2S = 1 
sets a threshold for the partonic reaction, where Q and 
yS denote the invariant mass of the lepton pair and the 
overall hadronic center-of-mass (cm.) energy, respec- 
tively, and x\ and X2 are the momentum fractions of the 
partons participating in the hard-scattering reaction. As 
z increases toward unity, little phase space for real-gluon 
radiation remains, since most of the initial partonic en- 
ergy is used to produce the virtual photon. The infrared 
cancellations between virtual and real-emission diagrams 
then leave behind large logarithmic higher-order correc- 
tions to the cross sections, the so-called threshold log- 
arithms. These logarithms become particularly impor- 
tant in the fixed-target regime, because here the ratio 
Q 2 /S is relatively large. It then becomes necessary to 
rcsum the large corrections to all orders in the strong 
coupling, a technique known as threshold resummation. 



QCD threshold resummation for the Drell-Yan process 
has been derived a long time ago 14). It turns out that 
the threshold logarithms lead to a strong increase of the 
cross section near threshold. Therefore, if threshold re- 
summation effects are included, it is possible that a much 
softer valence distribution of the pion is sufficient to de- 
scribe the experimental data. Indeed, as was observed in 
Ref. [6], the extracted v w already becomes softer when 
going from the lowest order to NLO, where the thresh- 
old logarithms first appear. In this Letter, we will ad- 
dress the impact of resummation effects on the pion's 
valence distribution. We will find that indeed a falloff 
jj" ~ (1 — x) 2 even at a relatively low resolution scale is 
well consistent with the Drell-Yan data. We note that the 
effects of resummation on parton distributions were also 
examined in the context of deep-inelastic lepton scatter- 
ing 0. 

We consider the inclusive cross section for the produc- 
tion of a /i + (i~ pair of invariant mass Q and rapidity 77 
in the process tt~ (Pi)A(P 2 ) — > /i + /i~A, where A denotes 
a nucleon or nuclear target and P\ and P2 are the four- 
momenta of the initial-state particles. According to the 
factorization theorem, the cross section is written as 



do 



00 



dW dx, A 

X 0J ab (xx,Xi,X 2 ,X2,Q/fJ,), (1) 

where a = 4ira 2 /9Q 2 S, with S = {P1+P2) 2 , and where 

7e ±7? (2) 



x l,2 



with t = Q 2 /S. At lowest order, one has X\^ — x i 2- 
The sum in Eq. ([1} runs over all partonic channels, with 
and fjf the corresponding parton distribution func- 
tions of the pion and the nucleus and u) a b the hard- 
scattering function. The latter can be computed in per- 
turbation theory as a scries in the strong coupling con- 
stant a s , starting from the lowest-order process qq — > 
7* — > /.i + fi~. The functions in ([I]) depend on the factor- 



ization and renormalization scales, which we choose to 
be equal here and collectively denote as fi. 

As discussed above, our goal is to resum large log- 
arithmic contributions to Lo qq that arise near partonic 
threshold z = 1, where z = Q 2 /s = t/x\X2, with the 
partonic cm. energy. Resummation may be achieved 
in Mellin moment space, where phase space integrals for 
multiple-soft-gluon emission decouple. For the rapidity- 
dependent cross section, it is convenient to also apply a 
Fourier transform in 77 16|, [17( . Under combined Fourier 
and Mellin transforms of the cross section, 



a(N, M) 
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(3) 



the convolution integrals in ([T]) decouple into ordinary 
products (l6l . 17 1 . Defining the moments of the parton 
distribution functions, 



A^EE / dxX»^f{x^% 



,2\ _ 



N-l . 



(4) 



and introducing the corresponding double transform of 
the partonic hard-scattering cross sections, 



w ab {N,M)= / dzz 



JV-l 



in 



df/e 



iM fj 



(5) 



where fj = r\ 
one finds 



a(N,M)=a ^f a ' N+ii 

a . 6 



hbi(xi/x2) is the partonic cm. rapidity, 



' /jf' W_i ~ r &ab (N, M) . (6) 



As was discussed in Refs. 17-19j], in the near-threshold 
limit z — > 1 or N — > 00, the dependence of u) a b=qq (N, M) 
on M becomes subleading and may be neglected. The 
resummed expression for Cb qq {N, M) then becomes iden- 
tical to that for the total (rapidity-integrated) Drell-Yan 
cross section and is given in the MS scheme by 



lni 



= Co 



Q 2 



,a s {^)\+2 / dC 



N-l 



- 1 



i-C 



r^ 2 Q 2 dkl A . 
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where A q (a s ) is a perturbative function, whose first two 
orders are sufficient for resummation to next-to-leading- 
logarithmic (NLL) order [l4j : 



A q (a s ) 



l A W 
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(8) 
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(9) 



Here C F = 4/3 and C A = 3. The first term in Eq. 
does not originate from soft-gluon emission but instead 
mostly contains hard virtual corrections. It is also a per- 
turbative series in a s , and we need only its first-order 
term: 

C g = ^C F (-4 + + I In Sp) + 0(a% (10) 

whose exponentiated form is given in Ref. [2lj. As was 
shown in Ref. [17| . rapidity dependence is slightly more 
faithfully reproduced if one shifts the Mellin moments to 
N ± iM/2 in u qq , which is a choice that we also adopt 
here. 



At NLL order, the expression in Eq. becomes [22 



In lu„ 



C q + 2h {1) {X)\nN + 2h (2) fx, 



(11) 



where N — Ne lB with the Euler constant jf, and 

A = b a s ((i 2 ) In A. (12) 

The functions hS 1 ' and collect all leading-logarithmic 
and NLL terms in the exponent, which are of the form 
a k s ln fc+1 N and a* ln fc A, respectively. They read 



2nb Q X 



[2A + (1-2A) ln(l-2A)], 



2) \ 2A + ln(l - 2A) 



2^ 2 b% 



^ (1 _ 2A) + ^_ M1 _ 2A)lll ^, 



(13) 
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where 



b = (UC A -2N f ) 

IZH 



(14) 



bi = 



24tt 2 



(17C A - 5C A N f - 3C F N f ) . (15) 



The resummed hadronic rapidity-dependent cross sec- 
tion is obtained by taking the inverse Mellin and Fourier 
transforms of Eq. ^ : 



da 



dQ 2 dri 



™e-^ [ C+i00 ^r-a(N,M). 
2?r JC-ioo 2tti 



(16) 

When performing the inverse Mellin transform, the pa- 
rameter C usually has to be chosen in such a way that all 
singularities of the integrand lie to the left of the integra- 
tion contour. The resummed cross section has a Landau 
singularity at A = 1/2 or N — exp(l/2a s 6o)j as a result 
of the divergence of the running coupling a s in Q for 
k± — > Aqcd- In the Mellin inversion, we adopt the min- 
imal prescription developed in Ref. 22[ to deal with the 
Landau pole, for which the contour is chosen to lie to 



the left of the Landau singularity. An alternative possi- 
bility is to perform the resummation directly in z-space 
[23| . We match the resummed cross section to the NLO 
one by subtracting the 0(a s ) expansion of the resummed 
expression and adding the full NLO cross section [17 1. 

The fixed-target pion Drell-Yan data 0, Q are in a 
kinematic regime where the partons' momentum frac- 
tions are relatively large, (x > 0.3), and hence the va- 
lence quark contributions strongly dominate. We can 
therefore only hope to determine the pion's valence dis- 
tribution v w = v% — d-1 — (F v = ul . Following 
the NLO Gluck-Reya-Schienbein (GRS) analysis [5|, we 
choose the initial scale Qo = 0.63 GeV for the evolution 
and parameterize the valence distribution function as 

xv*(x,Ql) = N v x a (l-xf(l + jx 5 ), (17) 

subject to the constraint J Q v*(x, Qfydx = 1. Since 
there is no sensitivity to the sea quark and gluon dis- 
tributions, we adopt them from the GRS analysis, ex- 
cept that we modify the overall normalization of the 
sea quark distribution so that the momentum sum rule 
^2i =q qg fa dxxfi(x) = 1 is maintained when we deter- 
mine the valence distribution. All distributions arc then 
evolved at NLO to the relevant factorization scale [i = Q. 

The free parameters in Eq. (|17p are determined by a fit 
to the pion Drell-Yan data from the Fermilab E615 ex- 
periment [H , applying threshold resummation as detailed 
in the previous section. The E615 data were obtained by 
using a 252 GeV ir~ beam on a tungsten target. We take 
into account the nuclear effects in this heavy target by 
using the nuclear parton distribution functions from Ref. 
24j . We use data points with lepton pair mass 4.03 GeV 
< Q < 8.53 GeV (between the J/ty and T resonances) 
and < xf < 0.8. Here, xp is the Feynman variable. In 
the near-threshold region, which is addressed by thresh- 
old resummation, we can use lowest-order kinematics to 
determine the relation between xp and the rapidity r\: 



Xp — x^ 



~o = V7sinh(77). 



(18) 



Since the E615 data have a nominal overall systematic 
error of 16%, we introduce a normalization factor K that 
multiplies the theoretical cross section. We find that the 
parameter 5 in (|17[) is not well-determined, and we hence 
fix it to 8 = 2, a value roughly preferred by the fit. In 
order to obtain a better picture of the physical content of 
our determined pion valence distribution, we perform fits 
for several different values of its total momentum fraction 
(xv K ) = Jq ot t (i, Qo)- Fixing {xif) makes one parame- 
ter in Eq. (|T7|) redundant, which we choose to be 7. We 
hence fit the remaining three free parameters a, j3 and K 
to the 70 data points using a \ 2 minimization procedure. 

The results are shown in Table HI for four different val- 
ues of the total valence quark momentum fraction 2(xv 7T ). 
One observes that fit 3 for which the valence carries 65% 
of the pion's momentum is preferred, with slightly higher 



or lower values also well acceptable. Most importantly, 
all fits show a clear preference for a falloff much softer 
than linear, with fits 2,3, and 4 having a value of (3 very 
close to 2. This is the central result of our work. The 
valence distribution xv 7 ' for our best fit 3 is shown in 
Fig. [TJ evolved to Q = 4 GeV. At this scale it behaves as 
(1 — x) 2 - 34 . Valence distributions obtained from previous 
NLO analyses [HIE|, which have a roughly linear behavior 
at high x, and from calculations using Dyson-Schwinger 
equations H, for which v* ~ (1 — x) 2A , are also shown. 
We note that for all our fits the factors K lie well within 
the normalization uncertainty of the data. 




FIG. 1: The pionic valence (v n ) distribution obtained from 
our fit 3 to the E615 Drell-Yan data at Q — 4 GeV, compared 
to the NLO parameterizations of [1] (Sutton-Martin-Roberts- 
Stirling) and [§| (GRS) and to the distribution obtained from 
Dyson-Schwinger equations Q. 

In Fig. [5] we compare the resummed Drell-Yan cross 
section obtained for fit 3 to some of the E615 data. We 
have chosen the factorization and renormalization scale 
fi = Q. As one can see from the figure and from Table 
HI the data are very well described. This also holds true 
for the CERN NA10 0] Drell-Yan data, which were not 
included in our fit, and to which we compare in Fig. [3] 
We also show the results obtained for our fit 3 when using 
only NLO (i.e., unresummed) partonic cross sections in 
the calculation. As seen in Fig. [31 these fall off too rapidly 
at large xp. The inset in Fig. [5] shows the variation of 
our NLO and resummed cross sections with /x, in terms of 
the quantity A M = [er(/i') — <r{n = Q)] /cr{n — Q), with 
p! = Q/2, 2Q. As one can see, the scale uncertainty 
is significantly reduced after resummation and becomes 
smaller than the experimental uncertainties, in particular 
in the region of high x p . This implies that our findings 
for the pion's valence distribution are stable with respect 
to the main theoretical uncertainty in the calculation. 

In conclusion, we have determined a new valence 
parton distribution function for the pion by reanalyz- 
ing pion-nucleon Drell-Yan data including thrcshold- 
resummed contributions to the cross section. The ob- 
tained valence distribution is much softer in the high- 



TABLE I: Results for our NLL threshold-resummed fits to the Fermilab E615 Drell-Yan data 



Fit 


2{xv n ) 


a 




7 


K 


X' (no. of points) 


i 


0.55 


0.15 ±0.04 


1.75 ± 0.04 


89.4 


0.999 ±0.011 


82.8 (70) 


2 


0.60 


0.44 ±0.07 


1.93 ±0.03 


25.5 


0.968 ±0.011 


80.9 (70) 


3 


0.65 


0.70 ±0.07 


2.03 ±0.06 


13.8 


0.919 ± 0.009 


80.1 (70) 


4 


0.7 


1.06 ±0.05 


2.12 ±0.06 


6.7 


0.868 ± 0.009 


81.0 (70) 
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FIG. 2: Comparison of our NLL-resummed Drell-Yan cross 
section based on fit 3 to some of the E615 Drell-Yan data. 
The inset in the upper figure shows the scale variation of the 
resummed and the NLO cross sections (see the text). 



x regime than that found in an NLO analysis, behav- 
ing roughly as (1 — x) 2 , in agreement with predictions 
from perturbative QCD and nonperturbative Dyson- 
Schwinger equation approaches. 

We are grateful to F. Yuan for useful communications. 
This work was supported by BMBF. M.A. was supported 
by a grant of the "Bayerische Eliteforderung." W.V.'s 
work has been supported by the U.S. Department of En- 
ergy (Contract No. DE-AC02-98CH10886) and by the 
Alexander von Humboldt Foundation. 



[1] J. S. Conway et al. [E615 Collaboration], Phys. Rev. D 
39, 92 (1989). 

[2] P. Bordalo et al. [NA10 Collaboration], Phys. Lett. B 
193, 373 (1987); B. Betev et al. [NA10 Collaboration], 
Z. Phys. C 28, 9 (1985). 

[3] For a review, see: W. J. Stirling and M. R. Whalley, J. 




-0.4 -0.2 0.0 0.2 0.4 0.6 



FIG. 3: Comparison of our resummed results for fit 3 to some 
of the CERN NA10 n~-W data 0] at pion beam energy 194 
GeV. We have used a normalization factor K — 1.045. 



Phys. G 19, Dl (1993). 
[4] P. J. Sutton et al., Phys. Rev. D 45, 2349 (1992). 
[5] M. Gliick, E. Reya and I. Schienbein, Eur. Phys. J. C 10, 

313 (1999). 

[6] K. Wijesooriya, P. E. Reimer and R. J. Holt, Phys. Rev. 

C 72, 065203 (2005). 
[7] G. R. Farrar and D. R. Jackson, Phys. Rev. Lett. 43, 246 

(1979); E. L. Berger and S. J. Brodsky, Phys. Rev. Lett. 

42, 940 (1979); S. J. Brodsky and F. Yuan, Phys. Rev. 

D 74, 094018 (2006); F. Yuan, Phys. Rev. D 69, 051501 

(2004) . 

[8] M. B. Hecht, C. D. Roberts and S. M. Schmidt, Phys. 

Rev. C 63, 025213 (2001). 
[9] T. Shigetani, K. Suzuki and H. Toki, Phys. Lett. B 308, 

383 (1993). 

[10] T. Frederico and G. A. Miller, Phys. Rev. D 50, 210 
(1994); A. Szczepaniak, C. R. Ji and S. R. Cotanch, Phys. 
Rev. D 49, 3466 (1994). 

[11] W. Melnitchouk, Eur. Phys. J. A 17, 223 (2003). 

[12] R. J. Holt and C. D. Roberts, Rev. Mod. Phys. 82, 2991 
(2010). 

[13] H. Shimizu et al, Phys. Rev. D 71, 114007 (2005). 

[14] G. F. Sterman, Nucl. Phys. B 281, 310 (1987); S. Catani 

and L. Trentadue, Nucl. Phys. B 327, 323 (1989). 
[15] G. Corcella and L. Magnea, Phys. Rev. D 72, 074017 

(2005) . 

[16] G. F. Sterman and W. Vogelsang, JHEP 0102, 016 
(2001). 

[17] A. Mukherjee and W. Vogelsang, Phys. Rev. D 73, 
074005 (2006). 

[18] E. Laenen and G. Sterman, FERMILAB Report No. 

CONF-92-359-T (unpublished). 
[19] P. Bolzoni, Phys. Lett. B 643, 325 (2006). 
[20] J. Kodaira and L. Trentadue, Phys. Lett. B 112, 66 

(1982). 



[21] T. O. Eynck, E. Laenen and L. Magnea, JHEP 0306, 
057 (2003). 

[22] S. Catani et at, Nucl. Phys. B 478, 273 (1996). 

[23] T. Becher, M. Neubert and G. Xu, JHEP 0807, 030 



(2008). 

[24] D. de Florian and R. Sassot, Phys. Rev. D 69, 074028 
(2004). 



